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Abstract: We study the long-Hved modes of the charge density and energy density 
' correlators in the strongly-coupled, finite density field theory dual to the AdS4, Reissner- 

D Nordstrom black hole. For small momenta q <^ fi, these correlators contain a pole due to 

I sound propagation, as well as a pole due to a long-lived, purely imaginary mode analogous 

to the /i = hydrodynamic charge diffusion mode. As the temperature is raised in the range 
^ T < fi, the sound attenuation shows no significant temperature dependence. When T > fi, 

it quickly approaches the /i = hydrodynamic result where it decreases like 1/T. It does 
\^ not share any of the temperature-dependent properties of the 'zero sound' of Landau Fermi 

liquids observed in the strongly-coupled D3/D7 field theory. For such small momenta, the 
^ energy density spectral function is dominated by the sound mode at all temperatures, 

whereas the charge density spectral function undergoes a crossover from being dominated 
by the sound mode at low temperatures to being dominated by the diffusion mode when 
J> T ~ fJ-'^/q- This crossover occurs due to the changing residue at each pole. We also compute 

^ the momentum dependence of these spectral functions and their corresponding long-lived 

^ poles at fixed, low temperatures T <C 
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1 Introduction 

The AdS/CFT, or gauge/gravity, correspondence [1-4] is a very useful tool which allows 
one to study the properties of certain strongly-coupled field theories via their dual grav- 
itational theories. In particular, the thermodynamic and near-equilibrium properties of 
such strongly-coupled field theories can be obtained relatively easily from their dual grav- 
itational descriptions. Initial studies of these properties concentrated on field theories at 
non-zero temperature T (most notably N" = 4 SU{Nc) supersymmetric Yang-Mills theory 
with Nc —7- oo) and were motivated by experimentally-observed properties of thermal field 
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theories [5]. For perturbations whose frequency w and momentum q are much less than 
T, these field theories were found to obey the laws of hydrodynamics and their transport 
coefficients such as shear viscosity, charge diffusion constant etc. were calculated (see [6-8] 
and subsequent work). More recently, there has been a lot of interest in studying field 
theories at zero temperature but with a non-zero density of a conserved global U{1) charge 
- these are analogues of strongly-coupled condensed matter systems with a non-zero density 
of particles (see [9-12] for some introductions to the field). 

There have been numerous studies of the excitations present when T = and uj,q <^ 
specific strongly-coupled field theories, and one common feature amongst many examples 
is a propagating longitudinal mode with the dispersion relation 

iTo<f + 0{(f), (1.1) 

and q=-. (1.2) 

This is a feature of probe brane theories in different dimensions and with different UV 
symmetries (where the conserved charge is a density of fundamental matter - at least in 
the case where the background is derived from string theory) [13-21] and also the 4D bulk 
Einstein-Maxwell theory with a cosmological constant (where the conserved charge is an 
R-charge density) [22]. 

In one of these probe brane theories (the D3/D7 theory in (3+1) dimensions), the be- 
haviour of this zero temperature sound mode was analysed as the temperature was increased 
from T = and it was found to behave similarly to the 'zero sound' mode due to the oscil- 
lation of the Fermi surface of a Landau Fermi liquid (LFL), despite the fact that some of 
its other properties are quite different from an LFL (e.g. its heat capacity is proportional 
to rather than T) [23]. The theory of an LFL is valid when T <C /i and w,g <C /i and 
it predicts three different regimes for the behaviour of the sound attenuation, as shown in 
figure 1 [24-28]. At T = 0, the sound mode has a non-zero attenuation V proportional to 
Lo^ I II. As the temperature is increased from zero (with w and ^ fixed), the attenuation 
of the sound mode remains approximately constant until T//i ~ w/^u (the 'collisionless 
quantum regime' denoted by A in figure 1). Above this, thermal excitations around the 
Fermi surface must be taken into account and collisions between these thermally-excited 
quasiparticles cause the sound attenuation to increase at a rate proportional to T^/// (the 
'collisionless thermal regime' B). When T is sufficiently high such that T//i ~ ^JujJJi., these 
thermal collisions become so frequent that the zero sound mode is no longer a long-lived 
mode. However, the thermal collisions support the hydrodynamic modes of sound and 
diffusion. This results in the 'collisionless/hydrodynamic crossover' and the sound atten- 
uation begins to decrease at a rate proportional to /i(a;/T)^ as the hydrodynamic sound 
mode stabilises. 

These features were reproduced precisely in the strongly-coupled D3/D7 field theory 
except that there was no sound propagation in the hydrodynamic regime - this can be 



uj = ±Vsq ■ 

where 

U}= -, 
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Figure 1: The sound attenuation F in an LFL as a function of temperature, at fixed uj and fi. A is 
tlie coUisionless quantum regime, B is tlie collisionless thermal regime and C is the hydrodynamic 
regime. 



explained by the fact that the fluctuations of the bulk metric, which generate the hydrody- 
namic sound mode of the dual field theory, were explicitly suppressed for consistency with 
the probe limit. The collisionless/hydrodynamic crossover is most clearly exhibited in the 
D3/D7 theory via the poles of the charge density Green's function in the complex frequency 
plane [23]. As the temperature is raised, the low temperature poles corresponding to sound 
propagation approach the imaginary axis and collide to form two purely imaginary poles, 
one of which becomes the hydrodynamic diffusion mode as the temperature is raised fur- 
ther. The crossover can also be seen by examining the charge density spectral function - 
the tall, narrow peak corresponding to the sound mode becomes smaller and moves towards 
the origin as the temperature is raised, eventually forming a diffusive peak around = 0. 
A natural question to ask is whether similar behaviour is observed in the low temperature 
sound modes of other holographic theories at large chemical potential and hence whether 
they can all be characterised in this way as LFL-like 'zero sound' modes or not. 

In this paper, we study the strongly-coupled field theory dual to the KN-AdS/^ black hole 
solution of 4D Einstein-Maxwell theory with a cosmological constant. Many properties of 
this theory have been investigated in recent years as it is relatively simple and yet has very 
interesting behaviour. When T = 0, the low energy behaviour of the field theory is governed 
by a CFTi dual to the AdS2 factor of the black hole's AdS2 x near-horizon geometry 
[29]. If one considers a probe Dirac action for fermions in this background, the field theory 
operators dual to these fermions exhibit Fermi surfaces of a non- Fermi liquid type [29-31]. 
If instead one considers a fermionic action which is the supersymmetric completion of the 
Einstein-Maxwell action, no such Fermi surfaces are observed [32-34]. 

The quasinormal modes of the bulk bosonic fields, which correspond to the poles of 
the retarded Green's functions of the field theory energy-momentum tensor T^'^ and U(l) 
current J^, have also been studied at zero temperature. At T = the transverse sector 
(i.e. transport perpendicular to the direction of the momentum flow) contains a branch 
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cut along the negative imaginary frequency axis, and no long-lived modes (by which we 
mean no modes satisfying a) — )• as g — ?■ 0). At very small temperatures T <C //, this 
branch cut becomes a series of poles, one of which has a dispersion relation of the form 
a) = —iVq^ + O (g'^) [35-37]. This is an analogue of the well-known hydrodynamic shear 
diffusive mode of the energy-momentum tensor [6]. At T = 0, the longitudinal sector (i.e. 
transport parallel to the direction of the momentum flow) also contains a branch cut along 
the negative imaginary frequency axis. In addition to this, it has two propagating modes 
with dispersion relations of the form (1.1) with Vg = {dP/ de)\T=o and Fq w fi'i](j/2{e + P), 
where r]o is the 'zero temperature viscosity' derived via the 'Kubo formula' 



lim — Im \G' 

w->0 oj '- 



R 



(1.3) 



T=0 



e and P are the field theory's energy density and pressure respectively, and the '~' signifies 
that these are equal to within about 10% [22, 38].^ A priori, this result is not obvious. It 
suggests that when T = 0, fj, acts as an 'effective hydrodynamic scale', at least as far as 
sound propagation is concerned. 

In this paper, we first study numerically the behaviour of the longitudinal poles for 
a fixed momentum ^ < 1 as the temperature is increased from T = to T S> /i. Of 
particular interest to us is the behaviour of the sound mode and whether there is a col- 
lisionless/hydrodynamic crossover as in the D3/D7 theory and Landau's theory of Fermi 
liquids. We find that the attenuation of the sound mode shows no significant temperature 
dependence over the range T < n where we may have expected to find LFL-like behaviour 
and above this it approaches the fi = hydrodynamic result [39, 40] where the attenua- 
tion decreases like T~^. Its speed is approximately l/\/2 at all temperatures. This is in 
complete contrast to the results for an LFL (as shown in figure 1). Such an outcome is not 
particularly surprising since all available evidence suggests that this field theory is not an 
LFL. However, this is also true of the D3/D7 theory and yet it possesses an LFL-like zero 
sound mode. Our results show that this kind of mode is not generic to all strongly-coupled 
field theories at large chemical potential which have a gravitational dual. 

The other long-lived, longitudinal mode of the theory corresponds to a purely imaginary 
pole that forms when the branch cut along the negative imaginary frequency axis dissolves 
at non-zero temperatures. This mode becomes more stable as the temperature increases 
and when T ^ its dispersion relation is that of the fi = hydrodynamic charge diffusion 
mode [40, 41]. 

In addition to the poles of the Green's functions, we also compute the spectral func- 
tions of energy density and charge density as the temperature is increased at fixed q < 1. 
The energy density spectral function is dominated by the peak corresponding to sound 
propagation at all temperatures. In contrast to this, the charge density spectral function 
undergoes a crossover from being dominated by the sound peak at low temperatures to 
being dominated by the diffusion peak at high temperatures T S> fJ-'^/q- Note that this 



^The sound attenuation To is only known numerically and thus an exact comparison cannot be made. 
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mechanism is quite different than in the D3/D7 theory where the sound poles cohide to 
form the diffusion pole. Here, the sound and diffusion poles coexist at all non-zero temper- 
atures (that we can access numerically) but their residues change and this results in the 
crossover. In the D3/D7 theory this crossover was reminiscent of that in an LFL, but here 
we know of no such comparison - in particular, the RN-74dS'4 crossover occurs outside of 
the 'quantum liquid' range T <^ fx where we may expect an LFL-like theory to apply. 

As well as the previous results regarding the temperature dependence of the poles and 
spectral functions at fixed momentum, we have also calculated their momentum dependence 
at various fixed temperatures. We find that the sound and purely imaginary modes exist at 
all non-zero temperatures that we can access numerically (T > 0.0219/i) with the dispersion 
relations 

u: = ^-iTiT/^i)f + Oif), (1.4) 

and 

u = -iD{T/fi)f + 0{f), (1.5) 

and we compute the functions D{T/ ji) and T{T/^) numerically. For this reason, and the 
fact that it becomes the ^ = hydrodynamic charge diffusion mode in the T ^ fi limit, 
we will refer to this purely imaginary mode as the diffusion mode. When q ^ 1, these 
modes no longer dominate the low-energy properties of the theory and we must consider 
additional poles of the Green's functions also. 

Finally, we investigate the properties of sound propagation in the theory, at some fixed 
momentum q, over the (T, fi) plane. As discussed above, it is known that when = 0, T is 
the hydrodynamic scale and sound will propagate provided its momentum satisfies q <^T. 
We also know that when T = 0, there is an effective hydrodynamic scale /x in that there is 
a long-lived sound mode provided that q <^ pL. We find that when both of these scales are 
non-zero, there will be a long-lived sound mode provided that any one of them is much 
larger than the momentum. In other words, sound will propagate for any value of /m/T", 
provided that one considers small enough momenta. 

The structure of the paper is as follows. In section 2 we present the RN-j4(i5'4 solution 
of the four dimensional Einstein-Maxwell theory with a cosmological constant, as well as 
the gauge-invariant fluctuations of the bulk fields that we use to compute the poles of the 
retarded Green's functions and the spectral functions. Section 3 contains the numerical 
results showing how the properties of the sound mode change as we increase the temperature 
at fixed q < 1. In particular, we demonstrate that these properties are significantly different 
from those of an LFL and of the D3/D7 theory. In section 4, we explore how the other 
long-lived mode of the system - the diffusion mode - behaves as T is increased for q <1. We 
also examine how the energy density and charge density spectral functions vary with T for 
q < 1. Section 5 contains results for the poles and spectral functions as a function of q at 
fixed temperature T < fi. We then study the existence of an effective hydrodynamic scale in 
section 6 by examining the properties of the sound mode as a function of both q/ fj, and q/T. 
This is followed by our conclusions and a discussion of the behaviour of the sound mode 
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compared to that of the D3/D7 theory in section 7. Appendix A contains the equations of 
motion and off-shell action for the aforementioned gauge-invariant fluctuations. 



2 The RN- 74^5*4 background and fluctuations 

2.1 The action and background solution 

The gravitational theory we study is the four dimensional Einstein-Maxwell theory with 
a cosmological constant, which has the action 

^ 1 " 



2nl 



[ (fxyf^ {R-2A- L^F^^Pi''') + 2 / d^x^/\h\ (/C + counterterms) 
Jm JdM 

(2.1) 

where A = — S/L^, h is the induced metric on the boundary of the spacetime, /C is the 
extrinsic curvature on this boundary and F^u is the field strength of a U{1) gauge field 
A^. This is a consistent truncation of = 11 supergravity [42, 43]. The counterterms are 
discussed in the appendix. 

This theory has a charged, asymptotically- AdS black hole solution with a planar horizon: 
the planar AdS^ Reissner- Nordstrom black hole (RN-A(iS'4) 

ly'^ f ( fY'\ ry'^ ry'^ T 2 

asA = yr^dt + — rdx + -TTjdy + „ „, . dr , 

r'^j^r) 

/M = l-a + <3^)(^)%Q=(^)\ (2.2) 



r 



where r is the bulk radial co-ordinate, rg is the position of the horizon and Q determines the 
C/(l) charge of the black hole. In the full D = 11 supergravity theory, it is the decoupling 
limit of the geometry created by a stack of M2-branes which are rotating (in a specific way) 
in the directions transverse to their world volume [43]. The bulk U{1) gauge field is dual 
to a U{1) R-current in the field theory on this worldvolume. 

This solution has one tunable dimensionless parameter Q which determines the ratio of 
the chemical potential in the field theory to the temperature 

It takes values between (the zero chemical potential limit) and \/3 (the zero temperature 
limit). The thermodynamics of the dual field theory are well-known [42]. In particular, we 
note that the entropy density has the unusual property of being non-zero when T = 0. 

2.2 Fluctuations around equilibrium 

We are interested in the response of the field theory to small perturbations around the 
equilibrium state - this is encoded holographically in the linear response of the black hole 
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to perturbations around the background solution (2.2): 

We use the rotational invariance in the (x, y)-plane to choose the momentum to flow only 
in the x-direction of the field theory, without loss of generality. We may then classify 
fluctuations according to their parity under y ^ —y. The fluctuations which are even under 
this operation {hxx, hyy, hrr, hu, hrt, hrx, hxt, o-r, o-t and a^) decouple from those which 
are odd {hyr, hyx, hyt and Uy) at linear order [41]. The indices are raised and lowered with 
the background metric. In the following we are interested only in the even fluctuations 
which we refer to as 'longitudinal' henceforth, as they encode the response of the fields 
parallel to the direction of momentum flow. The metric and gauge field fluctuations are 
coupled within this longitudinal sector which tells us that the retarded Green's functions 
of the longitudinal components of the field theory energy- momentum tensor T^'^ and U(l) 
conserved current are not independent. 

We are particularly interested in two properties of the retarded Green's functions Gq^^^. 
The first is the poles of the Green's functions in the complex frequency plane. These poles 
correspond to the field theory excitations - the real part of each pole is its propagating 
frequency and the imaginary part is its decay rate. We are primarily interested in the 
long-lived excitations - those with the smallest imaginary part. Note that if any excited 
bulk fields are coupled, their dual field theory operators share a common set of Green's 
function poles. 

The second object of interest to us is the matrix of spectral functions 

XAB q) ^ i {G^^os l) - GonO^ q)*) , (2.5) 

which tells us the rate of work done on the system by small external sources for Oa and 
Ob with frequency u (see, for example, [9]). Modes which couple strongly to external 
sources in this way are visible in the spectral functions as tall, narrow peaks centred on 
the propagating frequency and with a width proportional to their decay rate. Such a peak 
will be produced by a pole of the Green's function with small imaginary part provided that 
the residue of the Green's function is large enough at this pole and that there are no other 
poles near it in the complex frequency plane. Unlike the existence of a pole, the residue at 
a pole differs between the Green's functions of a set of coupled operators and hence despite 
the fact that they have a shared set of Green's function poles, the spectral functions of 
coupled operators can be very different. We are interested in the energy density spectral 
function Xt"T" = Xee and the charge density spectral function = XQQ^ which are real 
quantities, and tell us which modes couple strongly to external sources of energy density 
and charge density respectively. 
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2.3 Gauge-invariant variables and Ward identities 

The retarded Green's functions can be computed from the on-sheU action of the gravi- 
tational theory which generically has the form 



S. 



on-shell 



(l^lir, -w, -q)Gijdr4>j{r, uj, q) + 4>l{r, -OJ, -q)Cij(j)j{r, oj, q) 



(2.6) 

where (pj label the perturbations of bulk fields which are dual to field theory operators 
Oj. To obtain the retarded Green's function of two operators Gq^q^ one must 

find solutions to the bulk field equations which satisfy ingoing conditions at the black hole 
horizon and asymptote to 

(r — )• oo, — w, — g) — < ' and (r — ?■ oo, w, g) — < ' , (2.7) 

[O, I^A [O, J B 

near the boundary. Gq^q^ {^^q) is then given by the integrand of the on-shell action 
evaluated on these solutions [44, 45]. 

The longitudinal sector of our theory contains ten fields whose excitations are coupled. 
The equations of motion and on-shell action for the excited longitudinal fields can be 
simplified considerably by noting that the theory has a U(l) gauge symmetry which acts 
on the gauge field fluctuations as 

a^^a^- d^A, (2.8) 

and a diffeomorphism symmetry which acts as 

^^lu ^ hfj^i/ — V^i^^ — Vi/^^, 

^ «M - rv«^^ - AaV^e, 

to linear order, where V is the covariant derivative with respect to the background metric 
[46]. We can form two linearly-independent variables which are invariant under these 
transformations 

2r 



Zi{r,uj,q) = ioax{r,u,q) + qat{r,uj,q) ^ hyy, 



Z2(r,w,g) = -^hxt + -^hxx + -i^ritt H 2 — ""ra ^ + 



\ ^j{r) q j{r 

and to linear order in the fluctuations, we can write our theory in terms of these variables. 
In particular, this reduces the set of ten coupled equations to a set of two. It also allows 
us to write the on-shell action in the form 
didd'^q 



•Son-shell = / " 3" Zi{r,-u},-q)gijdrZj {r,u},q) + 4>i {r,-uj,-q)Cij4>j {r,uj,q) 

(2.11) 

where i = 1,2 and (pj denote the fundamental fluctuations {htt,hxx,hyy,htx,cLt,ax}- The 
Cij terms are analytic in w, q and hence contribute only contact terms to the retarded 
Green's functions (i.e. terms analytic in u},q). These Zi variables are a generalisation of 
those of [47] to non-zero chemical potential (and in 3+1, rather than 4+1, dimensions). 
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Written in this way we see explicitly that, neglecting contact terms, the on-shell action 
for a solution that has the form aj (r — >■ cx), zba;, ztg) — t- 1 (with all others fields zero in 
this limit) differs from that for a solution with (r — >• c«,iba;,ibg) — t- 1 (with all other 
fields zero in this limit) only by the factor g^/cj^ in the definition of the variables Zj, and 
similarly for the other fields. 

This property of bulk gauge-invariance thus generates a number of relationships between 
the retarded Green's functions of the corresponding field theory operators: 

2 

r<R — nR nR — ^ nR ^o\ 



nR —"^nR nR — nR nR — I t ^ \ nR T'i') 



2 

nR nR f^R _ ^ f~<R /o 1 4^ 

where the '. . .' represents other similar relations, and these equations should be understood 
to hold up to contact terms. These are precisely the Ward identities of the field theory. 

Thus not only do these gauge-invariant variables simplify the equations of motion for 
the bulk fiuctuations, they also explicitly encode the Ward identities of the field theory. 
The contribution of the contact terms to the diagonal retarded Green's functions is purely 
real, and thus they don't affect our results for the spectral functions Xee and XQQ- Note 
that the contact terms cannot be written in terms of these gauge-invariant variables - we 
believe that this is because the linear diffeomorphism transformations (2.9) are not the 
correct ones to apply to the action which is quadratic in fiuctuations. 

We note that these are not the only possible gauge-invariant choice of variables. Another 
choice is the Kodama-Ishibashi variables which involve radial derivatives of the bulk fields, 
and have the advantage that the two equations of motion in these variables decouple [48] . 

2.4 Equations of motion and on-shell action in dimensionless variables 

It is convenient to work with the dimensionless radial co-ordinate u = r/vQ. For T > 0, 
we use the gauge-invariant variables 



Zi (n, w, q) = UQx + qat - '^h^, 



Z2 {u,uj, q) = i/i 



2uqh^ + a;2/,^ - ffiu)hl + ff{u) [ 1 + ^ - ^ ] 



Uf'{u) LO 



2 



(2.15) 

where f{u) = 1 — {1 + Q'^)/v? + Q'^ /u^. The equations of motion and action in these 
variables are given in the appendix. To compute the poles and spectral functions, we used 
the numerical procedure described in [45]. 
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At r = 0, we use the Kodama-Ishibashi variables and follow the methods described in 
[22]. The equations of motion in these variables are given in appendix A of [22]. We could 
only obtain accurate numerics at T = above q > 0.1, and hence we only show T = 
results in this range. 

3 Temperature dependence of the sound mode 

Our primary motivation for studying this theory is that it supports stable, propagating 
excitations of energy and charge density at zero temperature and large chemical potential 
g <C 1. These sound modes at zero temperature have a dispersion relation of the form 
(1.1) where the speed is Vs = l/\/2 [22]. We want to know what effect the increase 
of temperature has upon this mode - in particular we are looking to see if it shares the 
characteristics of the 'zero sound' mode of a Landau Fermi liquid. This comparison can 
be made by studying the sound attenuation as a function of temperature for T ^ and 
looking for the three different regimes shown in figure 1. 

Note that when fj, = and uj,q <^T, there are sound modes with the dispersion relation 

At non-zero /x, we would expect to recover these results when fj, <^ uj,q <^ T, which is 
outside of the 'quantum liquid' regime T <C ^ where any LFL-like behaviour would be 
present. 

The temperature dependence of the real and imaginary parts of the sound mode are 
shown in figures 2 and 3 for various q < 1. Our finite temperature numerical results are 
shown along with T = numerical results (for q > 0.1 where we can obtain accurate 
results) and the /i = analytic result (3.1). 

The plots show that both the real and imaginary parts of the mode have a non-trivial 
temperature dependence. As the temperature is increased from zero, finite temperature 
corrections cause small changes to the real part of the sound mode whose sign depends upon 
the value of q. At sufficiently high temperature, T//x > 1, the real part quickly asymptotes 
to the ^ = hydrodynamic result (3.1). The imaginary part of the sound mode shows 
similar behaviour. This is slightly surprising - it indicates that the fi = result (3.1) is 
valid when q <^ fi <^ T. We will comment upon this again in section 6. 

To make an easier comparison with Landau Fermi liquid theory, we plot the temperature 
dependence of the imaginary part of the sound mode on a logarithmic scale in figure 4. 
These plots show only the region T < fi where we may expect such a theory to apply, and 
the imaginary part of the sound mode is normalised by Coq, its value at the lowest non-zero 
temperature we can access. 
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Sound mode, g=0.01 



Sound mode, ^=0.1 




0.208 



T//i _ 

Sound mode, q=03 





Sound mode, ^=0.5 



Figure 2: Variation of the real part of the sound mode as the temperature is increased. The crosses 
mark the T = numerical results, the dots are the numerical results for T > 0, and the solid lines 
are the n — analytic result (3.1). 



There is a stark contrast between these plots and the results expected for an LFL zero 
sound mode, shown in figure 1.^ Landau's theory predicts that as the temperature is 
increased at fixed q and fj,, the imaginary part of the zero sound mode should be approx- 
imately constant up until T//x ~ q/fi. Between T/fi ~ q/ and T/ji ^ ^ qj fx, it should 
increase like T^. Above T/^ ^ ^/qJJj- and below T/^ ~ 1, it should decrease like T^. 
None of these features are present in our results. The magnitude of the imaginary part of 
the sound mode in our theory shows no significant temperature dependence until T ~ //. 
Above this, it begins to approach the n = 0, to, q <^ T result (3.1) where it decreases as 
1/T. 



An explicit comparison between these RN-j4d54 results and the corresponding D3/D7 
results is shown in figure 5. This highlights the fact that the D3/D7 sound mode behaves 
like the LFL zero sound mode, whereas the JlN-AdS^ sound mode does not. 



^The magnitude of the frequency of the sound mode in the KN-AdSi theory is of the same order as its 
momentum for the results shown, and thus we compare to the LFL results with uj ^ q. 



(3 



Sound mode, g=0.01 



Sound mode, ^=0.1 




13 



1.0 1.5 

Sound mode, <7=0.3 




Sound mode, 9=0.5 




Figure 3: Variation of the imaginary part of the sound mode as the temperature is increased. The 
crosses marks the T = numerical results, the dots are the numerical results for T > 0, and the 
solid lines are the ji — Q analytic result (3.1). 



4 Further temperature-dependent properties of the theory when g < 1 

In addition to the sound modes at T = 0, there are other propagating modes lying deeper 
in the complex frequency plane as well as a branch cut along the negative imaginary 
frequency axis [22]. In this section, we study how this configuration changes as T//x is 
increased at a fixed momentum q < 1. In particular, we concentrate on the longest-lived 
purely imaginary mode - this exists at non-zero temperatures as the branch cut mentioned 
above dissolves into a series of poles when T 7^ 0. We note here that when |U = and 
u},q <^T, the longest-lived purely imaginary mode has the dispersion relation 

u; = -i^q' + 0{q^), (4.1) 

corresponding to hydrodynamic charge diffusion [41]. We expect to recover this behaviour 
at non-zero fi in the limit fj, <^ u,q <^ T. 

We also show in this section how the energy density and charge density spectral functions 
of the theory change with the temperature, and in particular how the residues of the long- 
lived modes play an important role in the transition from sound domination of the charge 
density spectral function to diffusion domination. 
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Figure 4: Variation of the imaginary part of the sound mode as the temperature is increased, in 
the regime T < /i. The dots are the numerical results for T > 0, and the two dashed lines on each 
plot denote T/ fi — q/ fi and T/fj, = \fqf\i as one moves to the right along the plot. 



4.1 Temperature dependence of the diffusion mode 



We begin by studying the longitudinal diffusion mode of the theory. At T = 0, the 
negative imaginary frequency axis is a branch cut [22]. At non-zero temperatures, this 
branch cut dissolves into a series of poles along the axis and generically these become less 
stable (they recede into the complex plane) as the temperature is increased. However, one 
of the modes is special in that it becomes more stable as the temperature is increased, 
and at very high temperatures it becomes the ;U = hydrodynamic charge diffusion mode 
(4.1). Figure 6 shows how the imaginary part of this mode changes with the temperature. 
Its real part is always zero. Unlike for the sound mode, this plot has the same shape for 
all values of 0.01 < q < 0.5 and so we show only one for brevity. The decay rate of this 
mode decreases monatonically as the temperature is increased, and is described well by 
the /i = result (4.1) when T > fi. Again, it is non-trivial that the // = result holds in 
the regime q <^ fj. <^ T. Note that there is no T = point on this plot because it does not 
make sense to ask where the pole is in that case - the whole negative imaginary frequency 
axis forms a branch cut. 
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Figure 5: A superposition of the plots of the temperature dependence of the normahsed imaginary 
part of the sound mode when = 0.2 for both the D3/D7 theory and the RN-A(iS'4 theory. 
Crosses denote the D3/D7 numerical results [23] and circles denote the RN-AdS'4 results. Moving 
from left to right, the dotted lines mark the transition points between the quantum and thermal 
collisionless regimes, and the thermal coUisionless regime and the hydrodynamic regime, in the 
D3/D7 theory. These occur when w ~ T and ui ~ /^J' respectively. There are no results for the 
D3/D7 sound mode in the hydrodynamic regime since the hydrodynamic sound mode is suppressed 
in the probe brane limit. We refer the reader to [23] for a more detailed discussion of these features. 



Diffusion mode, q-0.5 




Figure 6: Variation of the imaginary part of the longitudinal diffusion mode as the temperature 
is increased. The dots are the numerical results for T > 0, and the solid line is the ^ = analytic 
result (4.1). 



4.2 Movement of the poles in the complex frequency plane with temperature 

It is instructive to view the simultaneous movement of the Green's function poles de- 
scribed previously in the complex frequency plane as the temperature is increased - this is 
shown in figure 7. As described previously, the sound and diffusion modes both become 
more stable, approaching the real axis as the temperature is increased. Note that the sound 
mode is closer to the real axis than the diffusion mode for all values of the temperature 
and thus it is always the longest-lived mode of the theory. However this does not mean 
that it always dominates the low-energy properties of the theory, as we shall show in the 
following subsection. Finally, we note that both the sound and diffusion poles coexist for 
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Figure 7: Movement of the sound and diffusion poles in the complex frequency plane as the 
temperature is increased at fixed q — 0.5. An animated version of this figure is available at 
http: //www. physics . ox. ac.uk/users/Davison/RNAdS4animat ions 1 .html. 



all non-zero values of the temperature that we can access numerically (T > 0.02^). This is 
in contrast to the strongly-coupled D3/D7 field theory in which the low temperature sound 
poles collide to form the high temperature diffusion pole [23]. 

In addition to the sound and diffusion poles, there are 'secondary' modes corresponding 
to poles lying deeper in the complex frequency plane. These are much shorter-lived than 
the sound and diffusion poles when g < 1 and become less stable as the temperature is 
increased. They do not have a significant effect on the low energy properties of the theory 
when ^ < 1 and hence we do not show their temperature dependence for brevity. 
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4.3 Variation of the spectral functions with temperature 

Until now, we have focused exclusively on the positions of the poles of the retarded 
Green's functions in the complex frequency plane. While these are interesting, they do not 
tell us the full story of how charge and energy are transported in the theory. To investigate 
this, we now turn our attention to the spectral functions of the energy density and charge 
density. These determine the average work done on the system when an external source of 
some frequency is applied to either the energy density or the charge density respectively. 
Despite the fact that the retarded Green's functions of both of these operators have the 
same set of poles, their spectral functions are quite different as shown in figures 8 and 9. 

At very low temperatures, both spectral functions are dominated by the peak corre- 
sponding to sound propagation. As the temperature is increased, the spectral function 
of the energy density undergoes a fairly unremarkable change - the sound peak becomes 
narrower and taller (corresponding to a longer-lived excitation) but completely dominates 
at all temperatures. In contrast to this, the sound peak of the charge density spectral 
function becomes smaller (and narrower) as the temperature increases. At a sufficiently 
high temperature it becomes so small that the peak around a) = 0, corresponding to the 
high temperature diffusion mode, dominates the spectral function. At what temperature 
does this crossover occur? In figure 10, we show how the crossover value of /x/T (i.e. the 
value where the sound and diffusion peaks are of the same height) varies with q/ ji. There 
is a clear linear relationship, signifying that the crossover occurs when 



and that diffusion dominates when T ^ /i^/^- Note that since we are studying the range 
g < 1, this condition automatically implies that T ^ /x. This crossover is reminiscent of 
the /i — )• limit. In that limit, the fiuctuations of T^^'^ and decouple resulting in a 
hydrodynamic sound pole in the T^" correlators and a hydrodynamic diffusion pole in the 
jAt correlators. The charge density spectral function in this limit is shown in [49]. However, 
we cannot interpret the crossover shown above to be due to approaching this limit, since it 
corresponds to the limit /j, <^ q,T oi our results, whereas the regime we are studying here 
is g ^ ^ T. 

In the D3/D7 theory, the corresponding crossover occurred for Tcross. ~ \/qJ^ and was 
reminiscent of the collisionless/hydrodynamic crossover in an LFL [23], but we do not have 
a similar explanation here. In particular, the crossover observed here in the RN-^(iS'4 
theory occurs outside of the 'quantum liquid' regime T <C /i. 

5 Dispersion relations at fixed temperature T < fi 

In the previous sections, we established how an increase in temperature affects the sound 
and diffusion modes that exist at some fixed, low momentum q fi. We now turn our 
attention to studying the dispersion relations of these modes at a fixed temperature T and 
chemical potential This is done by fixing T/^ and varying q. 
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Figure 8; The energy density spectral function for q = 0.5 as the temperature is increased, 
in units of 2p?'r{)/K\. The peak due to sound propagation dominates at all temperatures. An 
animated version of this figure is available at http://www.physics.ox.ac.uk/users/Davison/ 
RNAdS4animationsl .html. 



5.1 The sound mode dispersion relation 

In [22], the dispersion relation at T = and q <^ 1 was found numerically to be of the 
form (1.1) with Vg = l/\/2 and Fq = 0.083, which is remarkably close to the dispersion 
relation expected from the 'zero temperature hydrodynamics' described in the introduction 
and in [22], which has Tq = 0.072 rather than 0.083. 

We found a dispersion relation of the form (1.4) to be valid for non-zero temperatures 
also. The quadratic coefficient of the attenuation F as a function of temperature is shown 
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Figure 9: The charge density spectral function for q = 0.5 as the temperature is increased, in 
units of 2ro/K|. There is a crossover between sound domination and diffusion domination at high 
temperature. An animated version of this figure is available at http://www.physics.ox.ac.uk/ 
users/Davisoii/RNAdS4cLriimationsl .html. 



in figure 11. These results were obtained by fitting over the range 0.01 < q < 0.5, and an 
example of the fit for T = 0.0219/x is shown in figure 12. At high temperatures T S> //, 
it agrees with the ^ = result (3.1), and at low temperatures T <^ ^ it approaches a 
similar value to the T = result of [22]. Although our results at very low T do not 
match smoothly onto the T = numerical result of [22], they differ only by around 10% 
and we believe that this is most likely caused by numerical inaccuracies, which grow as 
the temperature is lowered. The proximity of the numerical results to the prediction of 
T = hydrodynamics' is surprising - ultimately, an analytic calculation will be needed to 
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Figure 10: The dependence of the crossover value of n/T upon q/fi. The best fit straight-hne to 
these results has intercept ~ 0.003 and gradient « 0.34. 
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Figure 11: The temperature dependence of the quadratic term F in the imaginary part of the sound 
dispersion relation (1.4). Circles show our numerical results, the solid line shows the fj, = analytic 
result (3.1), the '+' shows the T = numerical result of [22] and the 'x' shows the prediction of 
T = hydrodynamics'. 
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Figure 12: The dispersion relation of the sound mode at T = 0.0219/i for 0.01 < q < 0.5. The 
circles show our numerical results and the solid line is the best fit w w q/V^ — i0.075q^ + 0{q'^). 
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Figure 13: The dispersion relation of the sound mode at two different temperatures: T 
(circles) and T = 0.159/i (crosses). The dashed line is the line Re (w) — q. 



determine whether these small discrepancies are due to inaccurate numerics, or whether 
this proximity is in fact a coincidence. The general trend of the results is clear however 

- as the temperature increases, the sound mode becomes more stable, as was observed in 
section 3. 

When q ^ 1, this series form of the dispersion relation is useless. Figure 13 shows 
the dispersion relation of the sound mode when q > fj, > T, at two different temperatures: 
T = and T = 0.159/i. The dispersion relations have the same shape at both temperatures 

- the real part asymptotes to Co = q, and the imaginary part tends to a constant, in the 
region q^ T. 

5.2 The diffusion mode dispersion relation 

Recall that at non-zero temperatures, the branch cut along the negative imaginary fre- 
quency axis becomes a series of poles and that the most stable of these becomes the = 
diffusion mode at high temperatures. Figure 14a shows the imaginary part of the disper- 
sion relation of this pole at two fixed, low temperatures T = 0.0219/i and T = 0.159/i (its 
real part is always zero). At both temperatures, the pole recedes quickly into the complex 
plane as the momentum is increased. Performing a polynomial fit to the imaginary part in 
the range 0.01 < q < 0.5 at the very low temperature T = 0.0219/i, we found a dispersion 
relation of the form (1.5) with D ~ 0.83. The fit is shown in figure 14b. This therefore is 
an analogue, at low temperatures T < q < fi, of the /i = 0, g <C T hydrodynamic charge 
diffusion mode. 

In fact, this quadratic form of the dispersion relation (1.5) is valid for all non-zero 
temperatures that we could access - the dependence of D upon T is shown in figure 15. 
We extracted D from a fit over the range 0.01 < q < 0.5. It decreases monatonically as 
the temperature is raised, in agreement with the results of section 4.1, and approaches the 
// = result (4.1) in the limit T //. Again, we note that this is despite the fact that we 
are studying the regime fj,^ q. 
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Figure 14: The dispersion relation of the diffusion mode at two different temperatures: T — 
0.0219^ (circles) and T = 0.159/i (crosses) with the polynomial best fit at T = 0.0219/i shown also 
(solid line). We cannot track the T — 0.0219// mode for as high momenta as the T = 0.159/i mode. 
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Figure 15: The temperature dependence of the quadratic coefficient D in the dispersion relation 
of the diffusion mode (1.5). The circles show our numerical results and the solid line is the analytic 
fi = result (4.1). 



We could not obtain the numerical accuracy required to access non-zero temperatures 
lower than T = 0.0219/i and hence we cannot say whether the mode exists with a quadratic 
dispersion relation for arbitrarily low non-zero temperatures. We emphasise again that this 
mode does not exist at T = itself (unlike the T = 'R-spin diffusion' mode of [16]), as 
there is a branch cut in that case. 

5.3 Dispersion relations of the secondary modes 

Figures 16 and 17 show the dispersion relations of the second stablest ('secondary') 
propagating and purely imaginary modes at two different temperatures. These differ qual- 
itatively from the sound and diffusion modes described above in that a) 7^ when ^ = 0. 
At high momenta q ^ ^ > T, the propagating secondary modes have the same form as the 
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Figure 16: The dispersion relation of the second stablest propagating mode at two different 
temperatures: T — {) (circles) and T = 0.159/j, (crosses). The dashed line is the line Re {uj) = q. 
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Figure 17: The dispersion relation of the second-stablest purely imaginary mode at two different 
temperatures: T = 0.0219/i (circles) and T = 0.159/i (crosses). 



sound modes - their real parts asymptote to Re (w) = q and their imaginary parts tend to 
a constant. The purely imaginary secondary mode recedes into the complex plane as the 
momentum is increased, and the rate at which this happens increases with the temperature. 

5.4 Movement of the poles in the complex frequency plane with momentum 

It is instructive to view the simultaneous movements of these poles in the complex 
frequency plane as the momentum is increased. This is shown in figure 18 for T = 0.159/x. 
We see that as q is increased, the purely imaginary modes both become less stable as 
described previously. This figure shows that the diffusion mode destabilises much quicker 
than the secondary imaginary mode. The propagating modes show a different behaviour - 
their speeds both increase but their imaginary parts move in opposite directions and begin 
to approach each other in the complex plane as q increases. They eventually cross, before 
moving off horizontally together along the relativistic trajectory Re (a;) = q. At these high 
values of q it is clear that our original separation of modes into the stablest (sound and 
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Figure 18: Movement of the six longest-lived modes in the complex frequency plane as a function 
of momentum, for fixed T — 0.159/i. The crosses denote the sound and diffusion modes, and the 
circles denote the secondary propagating and imaginary modes. An animated version of this figure 
is available at http://www.physics.ox.ac.uk/users/Davison/RNAdS4animations2.html. 



diffusion), second stablest etc. is of no value. 

5.5 Variation of the spectral functions with momentum 

Finally, we turn our attention to the spectral functions of the theory at low temperatures 
and as a function of the momentum q. These are shown in figures 19 and 20 for T = 0.159^. 
At g = 0.5, the lowest momentum shown, both spectral functions are completely dominated 
by the peak due to sound propagation. As the momentum is increased, this peak becomes 
smaller and wider and when q > fi, it no longer dominates the spectral function - a peak due 
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Figure 19: The energy density spectral function for T — 0.159/j, as the momentum is increased, 
in units of 2iJ?'rQ/n\. As the momentum is increased, the peak due to sound propagation becomes 
less dominant. An animated version of this figure is available at http://www.physics.ox.ac.uk/ 
users/Davisoii/RNAdS4cLriimations2 . html. 



to the secondary propagating mode also becomes visible. As the momentum is increased 
further, these two peaks merge into one peak which moves with speed Re (u)) = q and 
constant width. This is a direct reflection of the two corresponding Green's function poles 
approaching each other in the complex plane. At high momenta, the value of the spectral 
function is very small at low frequencies w < g and it only becomes significant when u) > q. 
Note that there is no significant difference between the charge density and energy density 
spectral functions in this regime. 
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Figure 20: The charge density spectral function for T = 0.159/z as the momentum is increased, 
in units of 2ro/Kl. As the momentum is increased, the peak due to sound propagation becomes 
less dominant. An animated version of this figure is available at http://www.physics.ox.ac.uk/ 
users/Davison/RNAdS4cLnimations2 . html. 



6 An effective hydrodynamic scale 



As discussed previously, when /x = there is a long-lived sound mode with momentum 
q provided that q <^ T, and when T = there is a long-lived sound mode provided that 
q <^ fi. In the first instance, this is the regime of applicability of hydrodynamics and the 
condition on the momentum is such that the perturbations occur over much larger distance 
scales than the mean free path between thermal collisions. 
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Figure 21: Contour plot showing |Ini((Ii) |/Re(w) for the sound mode as a function of q/fi and 
q/T and the best fit to these results: oq = 10.1 and ai = 8.3. Darker colours correspond to smaller 
values (i.e. more stable propagation) and the contours show the values 0.02, 0.04, 0.06, and 0.08. 



We have studied the behaviour of the sound mode when both T and are non-zero, to 
determine if there is some 'effective hydrodynamic scale' (or effective mean free path) which 
determines whether sound propagation is possible or not in this more general case. Figure 
21a is a contour plot showing the value of |Im (a)) |/Re (w) - which is the ratio of the decay 
rate to the propagating frequency - for the sound mode as a function of q/ fi and q/T. Darker 
colours correspond to smaller values (i.e. more stable propagation). There is a clear pattern 
in the plot - provided that one of q/T or q/^ is small enough, there is stable sound prop- 
agation. It suggests that there is an 'effective hydrodynamic scale' governing sound prop- 
agation which is qualitatively of the form E'eff. = T (ao + . . . + (fi/T)" -|- . . . + aifi/T) 
where a € (0, 1). This reduces to the correct form in the T = and /i = limits separately. 
A fit of the form q/(aoT + aifi) could not quantitatively reproduce the plot above - sug- 
gesting that this ansatz is an oversimplification (for example, it neglects almost all terms 
of the form (/i/T)" in the denominator as well as higher order terms in q) - but does give 
qualitatively the correct features for the sound propagation properties. Figure 21b shows 
the best fit to this form. 

An effective hydrodynamic scale of this form is also consistent with the fact that the 
region q <^ fi <^ T reproduces the ^ = 0, g <C T results, as we saw in section 3. 

7 Conclusions and discussion 

In summary, our main results are as follows: 
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i) When q <ti the long-hved modes of the charge density and energy density Green's 
functions are the sound and diffusion-hke modes with dispersion relations (1.4) and (1.5) 
respectively. 

ii) When q <^ fi and T fi, the attenuation of the sound mode shows no significant 
temperature dependence, unlike in the strongly-coupled D3/D7 field theory and in Landau's 
theory of Fermi liquids. When q <^ ^ and T ^ ^u, the sound and diffusion dispersion 
relations are well-approximated by the /U = results of [39, 41]. 

iii) When q <^ fi, the energy density spectral function is dominated by the sound peak 
at all temperatures, whereas the charge density spectral function undergoes a crossover 
from sound domination at low temperatures to diffusion domination at high temperatures. 
This crossover is due to the changing residue at each pole, and occurs at a temperature 
Across ~ /u^/g > n. 

iv) When q ^ n and T <C /U, the sound and diffusion modes no longer dominate the 
energy density and charge density spectral functions, and the effects of other modes become 
important. 

v) When both /x and T are non-zero, a long-lived sound mode will propagate provided 
that its momentum q is much less than either /i or T. 

Our results show that although many strongly-coupled field theories at large chemical 
potential (which have a gravitational dual) possess a T = sound mode, these are not all 
LFL-like 'zero sound' modes (by which we mean that they don't have the properties shown 
in figure 1) as was the case in the D3/D7 theory. However, there is still the possibility that 
there could be universal behaviour of the sound mode within subsets of strongly-coupled 
field theories with a gravitational dual. We note that the density-dependent physics in the 
D3/D7 and RN- 74^54 theories arise through different holographic mechanisms (see [50, 51] 
for further discussion of this). In the D3/D7 theory, the background metric is fixed and 
it is the gauge field action - the DBI action - which alters the equation of motion of the 
gauge field from the n = Maxwell equation (whose only long-lived mode is the high 
temperature charge diffusion mode). In contrast to this, the gauge field equation of motion 
in the RN-^d54 theory departs from the = Maxwell equation via couplings to the 
bulk metric fluctuations. The other major difference between these two field theories is the 
number of spatial dimensions, but we don't expect that this will have a significant effect 
on the acoustic properties of the theory (provided that the number of field theory spatial 
dimensions is greater than one). 

It would be very interesting to check whether the low temperature sound modes in other 
probe brane theories [17-21] share the LFL-like properties of the D3/D7 theory. This would 
help to establish whether it is the form of the DBI action, which implies that the non-zero 
density in the field theory is a density of fundamental matter (at least in the cases where 
the background geometry can be derived from string theory) rather than, for example, the 
R-charge density in the RN-^dS'4 theory, that generates these interesting properties or not. 
The effects of metric backreaction (i.e. coupling between the field theory's charge density 
and energy density) are also yet to be computed for these probe brane theories. These may 
complete the LFL-like picture of acoustic propagation (by reproducing the full LFL sound 
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attenuation curve - figure 1 - including the hydrodynamic regime), or they may result in 
deviations from it. We cannot gain any insight into this from our backreacted RN-^(iS'4 
results because of the difference in gauge field actions described previously. Although an 
expansion of the DBI action in powers of F^y yields the Maxwell action at lowest order, 
in field theory quantities this is an expansion in powers of /u/T which is the opposite limit 
from that in which any LFL-like properties of a theory would be exhibited. 

In addition to those mentioned above, there are numerous other field theories with a 
gravitational dual which may possess interesting sound modes when T <^ fi. Among the 
most interesting of these are dilatonic black holes [51-53] and geometries where the bulk 
charge density is sourced by fermions [54-58]. It would also be worthwhile to determine 
the acoustic properties of more general truncations of supergravity which admit more 
complicated solutions than RN-j4(i54 (for example, those of [43, 52, 59]), to determine 
whether the specific truncation chosen has a significant effect on these properties. 

We have seen that even this relatively simple holographic theory has many non-trivial 
features in its bosonic excitations. Among the most intriguing are the accuracy of the 
'zero temperature hydrodynamics', and the crossover temperature Tcross. ~ A*^/? between 
the charge density spectral function being dominated by the sound mode and the diffusion 
mode. It would be useful to have a clearer physical understanding of these properties, and 
also to determine if they are present in other field theories at non-zero chemical potential. 
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A Equations of motion and action 

This appendix contains the equations of motion and action used to compute our T > 
numerical results. Our T = results were obtained in the way described in [22]. 

A.l Equations of motion 

There are two linearly-independent, coupled equations of motion for Zi (u, Cj, q) and 
Z2 {u, Q, q) which can be written in the form 

Z'{{u) + AiZ[{u) + A2Z'2{u) + A^Zi{u) + A^Z^iu) = 0, 
Z'{{u) + BxZ'-^{u) + B^Z'^iu) + B^Zxiu) + B^Z2{u) = 0, 
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where we have suppressed the dependence of Z\^2 on u and g, and the coefficients are 



A. 



Ao = 



Bo 



Ba 



2 = -tq 



Q2[^2(^2__^(^)^-2)_4y.(^)] 4Q2^-2 



+ 



+ 



2(72 [u? 



fiu)f) + Q2^ 



n6 (cu2 - /(u)g2) 

2;721 



A4 = iq^ 
Bi 



(^2 _ /(u)g2) [uf'{u)(p - 4 (w2 - /(^/)g2)] ^ 
4Q2 + u^f{u) 



-iq 



[«/ W - 4 (a;2 - /(n)g2)] ' 
2Q2 [uf'{u)f-4{oj'-f{u)f)] 



(w2 _ /(ix)g2) 



3^ 



-16nV(«)'9 



S, = 



«5/(n) (cD2 _ ;(^)^2) [n/'(n)g2 - 4 (lD2 _ /(^)^2)] ^ 
+ 2f{u)^f [-4Q2g2 + i6txV + (/'(u) + uf'iu)) 
- f{u) -SQ^fuj'^ + 16nV + (2g2g2 _ 2^^4-2 ^ ^5 ^/^^^ -2^ 

+ 2u''f{u)oj^]] + - 4^;^] |^ 

SiQ'^q 



u^fiu) (^2 _ ;(,),-2) [,;.(,),-2 _4(c.2_/(,),-2)]|- / (-)'^ l^)^" " 1 

- /(«)g2 [-4Q2a;2 ^ ^_^/(^) (g2-2 _ 3^4-2 ^ ^5_^'(y)g2) + ^6_^//(^)-2] 



4QV(«) V + Q^^^ [uf'{u)q^ - 4a;2] 



nV(«)2[^/'W-4(a;2-/(«)g2)] ^ 
+ f{u)q^ [SQ^co^ + uf'iu) {u^fiu) + bu'f'iu) - Q^q^)] \ . 
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A.2 Action 

In these variables, the ofF-shell action to quadratic order in the fluctuations is of the form 



24 



du 



dujd^q 
7^ 



QuduZi{u, -w, -q)duZi{u, uj, q) + GuduZi^u, -u, -q)duZ2{u, u, q) 

+ G2iduZ2{u, -u, -q)duZi{u,u, q) + fe5„Z2(ti, -Q, -q)duZ2iu,uj,q) + 

(A.3) 
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where the coefficients are 





2u^ 


Qii = 








Ql2 = 








Q21 = 


-Ql2-, 



Q22 



- f{uW) I 
2u^f{u) 



2iquf{u) 
-A{u!^ - f{u)q'^)+uf'{u) 



(A.4) 



2 ' 



and the ' . . . ' represents terms with less than two u derivatives (which cannot generically be 
written in terms of these gauge-invariant variables). These coefficients, combined with the 
equations of motion listed previously, allow us to compute the Green's function's poles and 
spectral functions Xee, XQQ by following the method of [45]. Note that the counterterms 
in (2.1) (listed, for example, in [22]) do not affect these quantities. 

The numerical check described in [45] relies on the coefficients of the one-derivative terms 
in the action in addition to the two-derivative terms. Hence to obtain a numerical check 
in our gauge-invariant formalism, we added a boundary 'counterterm' (distinct from those 
mentioned previously) to the off-shell action of the form 



^C.t. 



24 



dujd^q d r 

du TT— (hj [u, 

{27rf du ^ 



-UJ, -q) Vfj {u, uj, q) (pj (n, w, q)] , 



(A.5) 



where "Dfj (u,oj,q) was chosen such that the full off-shell action could then be written 
in terms of the gauge-invariant variables Zi^2- This boundary term does not alter the 
equations of motion and as the coefficients Vfj are purely real, it does not have any effect 
upon the poles of the Green's functions or the diagonal spectral functions Xtt and XQQ- 
We do not list the coefficients here as they are very lengthy and do not directly affect the 
results presented. For our numerical computations, we found it more convenient to use the 
radial co-ordinate z = 1/u. 
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